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Abstract 

We prove that some non-self-adjoint differential operator admits factorization and apply 
this new representation of the operator to construct explicitly its domain. We also show 
that this operator is J-self-adjoint in some Krein space. 
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1 Introduction 

The time-evolution of a thin film of viscous fluid on the inner surface of a rotating cylinder can 
be approximated by the Cauchy problem for the periodic backward-forward heat equation 

Vt + l[y] = 0, y(0,x)=y , y(x,t) = y(x + 2n,t), x G ( 7r, 7r), t>0 (1.1) 

where 

l[y) = e(smxy x ) x + y x , e > 0. (1.2) 

This model was derived by Benilov, O'Brien and Sazonov [TJ [2] under assumption that the 
parameter e, related to the thickness of the film, is sufficiently small, i.e e << 1. 

It was first shown numerically by Benilov, O'Brien and Sazonov pQ that the spectrum of the 
operator L defined by the operation l[.] and periodic boundary conditions y{— tt) = y(ji) consists 
of pure imaginary eigenvalues only. This result was quite a surprising because it is well known 
that the Cauchy problem 11.11 is ill-posed at least for the class of finitely smooth functions. 
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The spectral properties of the operator L were studied rigorously by Davies [3] and by 
Chugunova, Pelinovsky |4j. Approaching the problem in two different ways they proved an- 
alytically that if the parameter e is within the interval |e| < 2 then the operator L admits 
closure in the Hilbert space L 2 [—n, ir] and being restricted to the orthogonal to a constant sub- 
space has compact inverse of Hilbert-Schmidt type and as a consequence its spectrum is discrete 
with the only accumulation point at infinity. 

The numerical conjecture that all eigenvalues are pure imaginary was recently proved by Weir 
[5]. The elegant proof based on the continuation of the eigenfunctions into Hardy space on the 
unit disk and explicit construction of the symmetric operator. 

Our goal in this paper is to find the factorization of the operator L and to construct its domain 
explicitly. As a consequence we prove that the operator L acting in physical space is unitary 
equivalent to the operator A introduced by Davies [3J acting in the Fourier space. We also prove 
that the non-self-adjoint differential operator L belongs to the class of J-self-adjoint operators 
in some Krein space. For basic facts related to Krein spaces see [5J. 

2 Factorization of the non- self- adjoint operator L. 

We denote by D(T) and 9t(T) the domain and the range of linear operator T respectively. The 
notation C 2 is used for the standard Lebesgue space of scalar functions defined on the interval 
(— 7T, 7r). From here on L is the indefinite convection-diffusion operator L : 



with the domain of all absolutely continuous 27r-periodic functions y(x) such that (Ly)(x) G C 2 . 
In addition, we define the operator S : 



(Ly)(x) = e ■ (sm(x)y' (x))' + y\x), C 2 i— > C 2 



C 2 ^C 2 , (Sy)(x) = y'{x), 
where y'(x) G C 2 , y(—n) = y(n), and the operator M : 

C 2 ^ £ 2 , (My)(x) = e ■ (sm(x)y(x))' + y(x) 
with the domain of all absolutely continuous functions y(x) such that (My)(x) G C 2 . 



Theorem 1 // the parameter e G (0, 2), then L = MS. 



Proof. Let us consider the operator A : 



C 2 I— > C 2 , (Ay)(x) = (sm(x)y(x))' 
with D (A) = {y(x) \ y(x), (Ay)(x) EC 2 }. Then a function y{x) can be written as 




(2.1) 
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If x > then 

4/2 



I 9(t)dt\ < ■ a(x) ■ x 1 

Jo sm(x) 



where a(x) = ( J* \9{x)\ 2 dx) 1 ^ 2 . Since the two summands in (12 .ip have different orders of growth 
as x — > this implies that if y(x) G C 2 then c = and 



sm(x) 

Moreover, 



y(x) = • / «(()<«. (2.2) 



A small modification of the same reasoning leads to the following estimation for every x G (— 7r, vt) 



2^|l/2 



with a(x) = | JJ^ |6*(x)| 2 (ix| 

Alternatively the same function y(x) can be written as 

y(x) = ■ (c - / 0(t)dt), 0(t) G £ 2 . (2.5) 

with the same 8(x) as in (12. ip . Representation (I2.5P yields the following relations 

y(x) = ^-- 1* 6(t)dt (2.6) 
sm(x) J x 

and 

ly(s)l< 1 ■ •/?(*) (2.7) 



sin(x) 



with/5(x) = |/J|0(a;)| 2 cfe| 1/2 . 

It follows from (E2D and (EE} that 



9{t)dt = 0. (2.8) 
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Starting from the point — ir one can also obtain that 



y(x) = ■ [ 9(t)dt, (2.9) 

(2 - 10) 



with 7(x) = (f* \6(x)\ 2 dx\ 1 ^ 2 and 

6(t)dt = 0. (2.11) 

Obviously, the natural domain of the operator B : 

C 2 i — ► C 2 , (By)(x) = sin(x) • y(x)', 

where y(x), (By)(x) G C 2 coincides with the domain of the operator A. And as a consequence 
we obtain that A* = —B. 

Now let us define the operator C : C 2 i— > C 2 , 

(Cy)(x) = -i ■ ( sin(x) • ?/(a;) - - cos(x)y(x)) = -i ■ ((By)(x) + - cos(x)y(x)), 

where y(x), (Cy)(x) G C 2 and D : 

(Dy)(x) = (| ■ cos(x) + 1) ■ y(x), C 2 ^ £ 2 . 

It follows directly from the relation between the operators A and B and inequalities (12.41) . (12.71 
(I2.10p that C is a self-adjoint operator. 

We now restrict the parameter e to the interval < e < 2. This implies that D has the 
bounded inverse, and hence that D^ 1 ^ 2 ■ C ■ D^ 1 ^ 2 is a self-adjoint operator. Therefore the 
operator defined as 

(i-e. D- 1 ' 2 ■ C ■ D- 1 ' 2 + /) 
also has the bounded inverse, and the same is true for the operator 

M = (i-e-C + D) = e-A + I. 

Finally, let us show that the subspace {const} 1 - n C 2 is invariant under M and M~ l . Indeed, if 
y(x) belongs to the domain of M, then y(x) has representation (12.21) . and hence 

(My)(x) =y{x)+e-9{x). 

It follows by (EHD and fl2TTTD that 

f 7T /»7T 

(My)(x)dx = / y(x)dx, 



so (My)(x) G {const} 1 - fl £ 2 if and only if y(x) G {const} -1 D £ 2 . 

Our main goal is to show that L = MS. It is easy to check that a function y(t) belongs to 
the domain (within C 2 ) of L only if 

y'(x) = (tan(|z|/2))- 1/e sin- 1 (x) • (c+ / (tan(|t|/2)) 1/e 0(t)rft), (2.12) 



where <f>(t) G C 2 . Let x > 0. Then 



(Un(t/2))^cf>(t)dt\ < ( / t^dtf /2) -(l |0(i)(^M)Vf d t)a/2). 
'o Jo Jot 

The latter estimation yields 

| f (t&n(t/2)) l/t <p(t)dt\ < x 1/2+1/t ■ a(x), 
Jo 

where — > if x — > 0. Then 

|(tan(x/2))- 1 / e sin- 1 (x) • ( / (tan(t/2)) 1/e 0(t)rft)| < x~ 1/2 ■ (3(x), 

Jo 

where (3{x) — > if x — > 0. If c = 0, we obtain that 



— 2/(0)| = | / (tan(r/2))- 1/e sin- 1 (r) • ( / (tan(t/2)) 1/e 0(t)dt)dr| < 7 (x)x 1/2 , 
Jo Jo 

where 7(2) — > if a; — > 0. Thus, if c = 0, is continuous at zero. At the same time if c 7^ 0, 
y(x) contains an additional summand of the order x~ l l e that is out of C 2 for e < 2. Thus, c = 
in f[2TT2"j) and so 

y'(x) = (t&n(\x\/2))- lle sin- 1 (x)- I (tan(|t|/2)) 1/e 0(t)rft . (2.13) 

Jo 

This representation was derived under the hypothesis that x > 0, but it is clear that it is valid 
for every x G (— 7T, 7r). 

Now let us assume that there is y(x) such that 

y(x)e®(L) but G'S(MS). (2.14) 

Then there are two options. 

• (Ly)(x) g m(MS) = {const} 1 - n C 2 ; 

• (Ly)(x) G m(MS). 

The first option means that 9\{L) = C 2 and without loss of generality one can assume that 
(Ly)(x) = 1. Then by (12131) 

y\x) = {t&n(\x\/2))- l / e s\n-\x) ■ AtanQ^)) 1 /^ . 







Thus, y'(x) > for every x G (— 7r, 7r). The latter is impossible for absolutely continuous 
27r-periodic function. It is a contradiction! 

Now let us consider the second option. If y(x) satisfies f)2.14p and (Ly)(x) G 9\(MS), then 
there is z(x) G D(MS) such that (Lz)(x) = (Ly)(x). The latter yields (L(z — y))(x) = 0. Thus, 
by virtue of (12.131) . y(x) = z(x) + const that is impossible thanks to (I2.14p . □ 



5 



Corollary 1 L is a closed operator with the non-empty resolvent set and its resolvent has Sturm- 
Liouville property. 

Let us define by L the operator that represents the restriction of L on the set of all 27r-periodic 
smooth functions and by L its closure. 

Theorem 2 L = Lq. The operator L is J -self- adjoint in the Krein space with indefinite metric 
J defined as J(f(x)) = /(-7T — x). 

Proof. It is evident that L C L. In the paper of Davies [5] was shown that L Q has the non- 
empty resolvent set. If L ^ L then every A G p(L ) belongs to cr p (L) but it is impossible due 
to Corollary [TJ 

The adjoint operator is defined by the operation 

(L*y)(x) = e • (sm(x)y'(x))' - y'{x) 

on the same domain as the operator L. The last statement of the theorem follows immediately 
from the equality L = JL* J . □ 

Remark 1 The J -self- adjoint operator L being restricted to the subspace orthogonal to a con- 
stant has a compact inverse that implies that there exist two, non-positive and non-negative with 
respect to the metric J, maximal invariant under the operator L subspaces. 

3 Domain of the operator L 

As it was proved in the previous section, the operator M : 

C 2 i— > C 2 , (My)(x) = e ■ (sm(x)y(x))' + y(x) 

with the domain of all absolutely continuous functions y{x) such that (My)(x) G C 2 has bounded 
inverse. Hence, there is a constant p\ such that 

||e • (sm(x)y(x))' + y(x) \\& > p x ■ \y{x) \ h i (3.1) 

for every y{x) G D(M). Next, the operator L is closed, so its domain is closed with respect to 
the norm of the graphic of L, i.e. for the norm 

l|y(z)ll, := {Ily(*)lli> + \\e-(sm(x)y'(x)y + y'(x)\\ 2 L2 y/ 2 . (3.2) 

Our aim is to show that for y(x) G Q(L) fl {const} 1 - the norm (13.21) is equivalent to the following 
norm 

\\y(x)\\ m := {|b'(x)||| 2 + || sin(x) ■ y\x)f L , + ||(sin(x) • y'(x)Y\\ 2 L ^ 2 . (3.3) 
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Indeed, for every y(x) G fl {const} 1 - we have 

\\y(x)\\ g < {\\y'(x)\\l 2 + 2||e ■ (sin(x)y'(x))'|| i2 + 2 • ||y'(x)||| 2 } 1/2 < p 2 ■ \\y(x) 
where p 2 = max{3, 2e}. From the other hand, taking into account (13.11) . we have 

\\y(x)\\ g > \\e- (sm(x)y'(x))' + y'(x)\\ L 2 > 



Pi 
2 



3^3 ■ \\y'(x)\\ L 2 +p 3 ■ \\e ■ (sm(x)y'(x))' + y'(x)\\ L 2 > 
P3 ■ (\\y'(x)\\ L 2 + || sm(x)y'(x)\\ L 2 + ||e ■ (sin(x) ■ y'(x))'\\ L 2), 
where p% = mm{^, |}. 

It follows from above that the domain of L is the linear sub-manifold H of the Sobolev space 



®(L) = H:feH 1 ( 



-7T, 7T 



/(7r) = /(-7T,) sm(x)f eH\-n,n) 



and is a Hilbert space with the norm defined as: 



||/'|& + ||sin(s)/'(s)|| 



H 1 



Remark 2 As a consequence of Theorem^ the tridiagonal matrix operator A in Fourier space 
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.4 



admits the same type of the factorization as in Theorem [I] 


e 

.4 BC= -3/2e 1 3/2e 

-2e 1 
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e/2 


— e 


1 





-3/2 









10 

2 

3 

4 



(3.4) 



(3.5) 



there B has bounded inverse and C has inverse of the Hilbert- Schmidt type. The operator A is 
J -self- adjoint with J = diag(— 1, 1, — 1, 1, ...). Thus, the natural domain ^£>{A) is the Hilbert 
space of all number sequences {/ ra }i° with the norm 



\\{fn}T 



{ $> 2 (l^| 2 + K n + - (" - l)/-i| 2 )} 



2 U l/2 



n=l 
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In the latter formula we put for definiteness /o = 0. Note also that there exist two, non-positive 
and non-negative with respect to the indefinite metric defined by J , maximal invariant under 
the operator A, subspaces. 
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